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On the Computation of the Far Field of 
Open Cassegrain Antennas 

By H. ZUCKER 

(Manuscript received February 2, 1970) 

The far field radiation patterns of open cassegrain antennas have 
been previously computed by a double integration 1 . Recently one 
integration has been eliminated by using a modified stationary phase 
approximation 2 . The rather remarkable accuracy thus obtained moti- 
vated an investigation to estimate the error in this approximation. 

In this note we present an alternate formulation for the computation 
of the far field and provide an additional estimate for the error in the 
modified stationary phase approximation. It is also shown that for 
certain illuminations, there exists a closed form solution to the azimuthal 
integration. For illuminations derived from a linear combination of 
TE°„ and TMJ„ waveguide modes, the integrals can be expressed in 
terms of convergent series of Bessel functions. 

Based on the aperture field method, the far field electric E f of an open 
cassegrain antenna is related to the reflected field at the aperture, E r , 
by the following double integral 1 



Bi _ ,- " p( - jtBJ jr'"jf &(«,») 



• exp [jk sin 0„(.r p cos 4> a + y p sin <£ a )]r 2 sin dd d<f> (1) 

where R a , U , <t> a are the far field spherical observation coordinates; 
X is the free space wavelength and k, = 2ir/A is the propagation constant; 
and m is the illumination angle. Referring to Fig. 1, the rectangular 
x p and y„ coordinates are related to the spherical 0, <f> coordinates by 

x„ = r(cos 0„ sin 8 cos cf> + sin 8 cos 6), (2) 

y„ = rsin 0sin0, (3) 

r = 2f/(a - b cos 0) , (4) 
where 

a — 1 + cos 6 cos O , (5) 
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Fig. 1 — Open cassegrain antenna. 
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b = sin e sin O , (6) 
and 

a 2 - b 2 = c 2 , (7) 

c = cos 6 + cos d . (8) 

d is the offset angle, and / is the focal length of the paraboloid. 

For combined TE° n — TMJ„ mode excitation and for y polarization 
of the horn feed the reflected electric field at the aperture, (E ru ), can 
be written 

(1 + cos d )(E ru ) = l ip sin J cEM cos + (6 + a cos <f>)E c (^ 



- u[(a 



cos - b)E e (0) cos <j> + eB.fS) sin2 (9) 



where E c (0) and E c (ir/c) are the fields radiated by the subreflector in 
the principal planes <f> = and = r/2 at a distance 2//(l + cos O ) 
from the focal point. 

The <£ integration in equation (1) was recently approximated by a 
modified stationary phase approximation 2 . An examination of this 
approximation reveals that the resulting phase factors and the arguments 
of the Bessel functions have a simple geometric interpretation. Specifi- 
cally the phase factors are directly related to the location in the x p , y v 
plane of the centers of the 6 = constant circles (see Fig. 2 in Ref. 2) 
and the arguments of the Bessel functions are related to the radii of 
these circles. 

This suggests, for the <f> integration, the following coordinate trans- 
formation 

x„ = r c cos 0! + x c , (10) 

y p = r e sin0, (11) 

with 

r. - 2^L» , (12) 

x. = S^aA. (is) 

This transformation gives in the x v , y v plane, for the 6 = constant 
coordinates the same set of circles as obtained from equations (2) and 
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(3), and for the 0i = constant coordinates a set of hyperbolas with the 
equation 

x 2 p - yl - 2x p y P ctn X + 4/(z, - y„ ctn 0.) ctn 6 n = (2/) 2 . (14) 

The sets of curves which correspond to the transformation (10) and 
(11) are shown in Fig. 2. 

In the 0, 0i coordinate system the rectangular x, y components of 
the far field have a much simpler representation. In particular for y 
polarization these components are 

(M,,)- ^^™ £'($,) si* ede (is) 

where the subscript y is to designate the y polarization, and the x, y 
components of (E u ) are 

«• = -(l+tos^c exp <*■ sin '■ cos *■) 

• f E c (0)(a cos 0! + b) cos 0! + #>(|)« « in2 0i 

exp [jfer c sin O cos (& - 0,, )] ^ ^ 

a + b cos 0! 

and 

. _ if exp Q'fc.r c sin 0„ cos 0„) 
( J ' " (1 + cos 0„)c 

• / E c (0)(a cos 0! + 6) - Er\V)c cos 0, sin fa 

exp [jfoy sin 0„ cos (0! - 0„)] , ^ 

a + 6 cos 0! 

In deriving equations (15), (16) and (17), the Jacobian, J, of the trans- 
formation (10) and (11) was used, which is 

(2i) 2 
J = Ht~ (a + b cos 0,) sin 0. (18) 

It is readily shown that by using a two term expansion of 

[l/(o + 6cos0,)] W (l/o)(l - 6/OCOS0O, 

the integrals (16) and (17) give the same values as those obtained by 
the modified stationary phase approximation but without the term 
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Fig. 2 — Transformation of (0, <£i) coordinates in the (x p , y p ) plane. 



(sin 6 sin 6 cos 0„) 2 in the denominator. This shows that the modified 
stationary phase approximation is at least of order (b/a)~, for all ob- 
servation angles and with the additional property that for large angles 
off axis it is of order [(b/a)'/'2kf sin 6.,]. 

There are two cases for which the integration of either equations 
(16) or (17) can be performed in closed form. 



(0 



aE e (0) = cE f \£\ 



where a and c are given by equations (5) and (8). 

Such illuminations can be approximately realized by combined 
TE°i — TM?! mode excitations. For this case the integration of equation 
(16) yields 

(w \ q M ex p (jkx f sin fl„ cos 0„) , 

{EX - -S* J (1 + cos do)c E c (0)Jo(x) (19) 
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where J is a Bessel function of order zero with the argument 

x = kr c sin 9 B . (20) 

The evaluation of the corresponding expression based on the modified 
stationary phase approximation yields 

2/2 



.„ . 87T 2 / 2 exp (jx c sin d a cos <f>„) 

C v)v ~ (1 + cos d )c 



•^0)L(,)+^r^-^^cos2^ 



(21) 



The error in the approximation is the second term in equation (21). 
That term is zero in the planes <f> a = tt/4 and <£„ = 3/4tt. In other planes, 
the error is of order (b/d) 2 /x, and is in agreement with the error estimate 
given above. 

The dependence of the second term in equation (21) on is shown 
in Fig. 3 for 6 = 55°, / = 80X and d a = 2.5°. These parameters are the 
same as those used in the computations in Ref. 2, Fig. 5. A comparison 
of the two figures shows good agreement, as would be expected, since 
the condition for the illumination was approximately satisfied. 

iii) cEM = aE^)- 

These illuminations may also be approximately realized. For this con- 
dition the cross-polarized component has a closed form solution given by 

(EX = ^ 2 / 2 ^P(fe s ^^ C0S ^) Ec (0) (-) JM sin *. . (22) 
v v (1 + cos d )C w 

An evaluation of the corresponding expression based on the stationary 
phase approximation gives* 



,„ . . 8tt 2 / 2 exp (jkx c sin 6 a cos &,) 

(£,„), -J (1 + cog e ^ c 



*»® 



J,(x) sin O + -2 ,2 n 2 , ~~ sin 3<^ a 

v a — o cos <b a x 



(23) 



The error in the modified approximation is the second term and is 
of the same order as the error term in equation (21). 

For illuminations derived from a linear combination of TEj n and 
TMJ, waveguide modes, the integrals (15) and (16) can be expressed 
in convergent series of Bessel function, by using the following Fourier 

* Equation (21) in Ref. 2 has a misprint and Bhould be divided by 2. 



OPEN CASSEGRAIN ANTENNAS 



1601 



1.25 



O 0.75 - 



S 0.50 - 



0.25 




16 20 24 

6 IN DEGREES 



Fig. 3 — Error introduced by the modified stationary phase approximation. 
Angle off axis, a = 2.5 degrees. <f> a = (TI plane); boo^ = v /2 (E plane). 



series expansions 

00 

exp [jx cos (4> - <£„)] = J (.r) + 2 X) f «f»C*0 c °s »(* _ 0«) 

n=l 

and 



1 



a + & cos 



= « + X) c n c °s ncf) 



with 



2 r 

!« = ~ / — 

7T .'„ fl 



COS 7fc£ 



rf0. 



-+- 6 cos <f> 
The integral (26) is tabulated 3 and by using equation (7) gives 

c ' \a + c. 



(24) 
(25) 
(26) 

(27) 



With these expansions, after some mathematical manipulations the 
following expressions are obtained for the integrals 



~ S "i e ?! J l c !u^ C T - {c[EM + E c \ :, 



(1 + cos d,y)(a + c)c 
+ [aEM ~ oE e (^j [i^^V.Cr) cos*,, 



/„(*) 



2c 



a — c 



„ c E(-n^) ^Wcos( 



n + 2)</>„J} 



(28) 
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and 
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fjp x _ 8tt/ 2 cx p (jx r sin fl„ co s 6.) J ., 
(*.)- - (1 + cos «„)(a + cfc V" 



(a + 2c)E c (0) + cEA- 



-./,(.,-) sin <£„ + 2c[c# e (0) - aE e (|) 

• £ (-i)"(^-:) n/2 - / n +2 W sin (n + 2)<f> 

n = () \t# ~t~ C/ 



From eguations (r>) and (7) 



5ju 



a — c 

g + c 



= tan 5/2 tan 0„/2. 



(29) 



(30) 



For most antenna designs, equation (30) is much smaller than one, 
the series in equations (28) and (30) are, therefore, rapidly converging, 
and can in principle be evaluated to any desired accuracy. 

It is noted that in the principal planes, these series can be related 
to Lommel functions of two variables . 

The stimulating discussions with N. Amitay and E. R. Nagelberg 
are gratefully acknowledged by the author. 
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Holographic Thin Film Couplers 
By H. KOGELNIK and T. P. SOSNOWSKI 

(Manuscript received June 17, 1970) 

Recently P. K. Tien and his co-workers have described a prism 
coupler as a convenient means to feed light into a single mode of a 
guiding thin optical film. 1 Distributed couplers of this kind are of great 
interest for integrated optical devices. In this brief we describe thin 



